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ABSTRACT 


This report provides a brief description and running instruc- 
tions for the one-dimensional Monte Carlo code ZTRAN. The program 
is used to calculate the transport of electrons and photons in 


heterogeneous multi-layer media. 


1. Introduction 

This report pertains to the Monte Carlo program code ZTRAN which has been 
developed as a tool for calculating the transport of electrons and photons 
through multi-layer media. The code is a generalization of an earlier Monte 
Carlo code, ETRAN,!-9 which applies to a homogeneous medium. In ZTRAN, the 
medium is assumed to consist of several adjacent plane-parallel layers, each 
of which can have a different composition. The layers are assumed to be 
unbounded laterally, and the treatment is one-dimensional, including only a 
Single spatial variable along the direction of stratification of the medium. 
The output of the ZTRAN code consists of information about the spatial distri- 
bution of energy deposition and the electron and photon fluxes (differential 
in energy) within the multi-layer medium, and the reflection and transmission 


of these radiations from the medium. 


Some results with ZTRAN have already been published. !9-1!2 The purpose of 
this report is to provide information which will facilitate the use of the 
code by others. Brief descriptions are given of the Monte Carlo model, and of 
the input and output of ZTRAN. Illustrative results are shown pertaining to 
energy deposition. Appendix A contains running instructions, in the form of a 
detailed listing of the input variables including their format and meaning. 
The purpose of the various subroutines in ZTRAN is also briefly indicated. 
Appendix B presents a listing of the input and output for a typical ZTRAN 


run. 


2. Monte Carlo Model 

The ZTRAN program is based on the same Monte Carlo model, and uses the 
same cross sections, as ETRAN. In fact, it shares with ETRAN the data base of 
pre-processed cross sections provided by the DATAPAC code.2 Successive photon 
interactions are sampled individually in direct analogy to the physical 
processes (Compton scattering, photoelectric absorption, and pair production). 
The much more numerous elastic and inelastic interactions between electrons 
and atoms are treated according to a “path-segment" model. The electron 
trajectories are divided into many segments in each of which numerous interac- 
tions occur. The net angular deflection from the combined effect of the 
elastic and inelastic collisions in a single segment is sampled from the 
Goudsmit-Saunderson multiple-scattering distribution, using — as the under- 
lying elastic scattering cross section — the Mott cross section modified for 
screening according to Moli@re. The total energy loss in a segment is 
sampled from the Landau straggling distribution, and from the appropriate 


bremsstrahlung cross sections. 


The Monte Carlo model takes into account all secondary radiations, 
including knock-on electrons from electron-impact ionization events, 
bremsstrahlung and characteristic x rays, Compton electrons, photoelectrons, 
electron-positron pairs, and annihilation radiation. The calculation follows 
all generations of electrons and photons in the target with energies above 
chosen cut-off values. A weighting scheme allows for the sampling of 
bremsstrahlung photon histories in excess of the natural production rates so 
that statistical fluctuations in the bremsstrahlung results can be reduced 
without an increase in the number of time-consuming electron histories. The 
statistical accuracy of the photon results is also improved by scoring for 
each collision point in the photon history the probability that the photon 
crosses the boundaries of interest without further interaction (collision- 
density method). 

A number of options can be selected which allow for variations of the 
basic Monte Carlo model. These are included to permit model studies involving 
various approximations, and are not necessarily recommended for production 
runs. These variations include, for example, treating electron energy losses 
in the continuous-slowing-down approximation, not following energetic knock-on 
electrons (delta-rays), and not following secondary photons (bremsstrahlung, 
characteristic x rays, annihilation radiation). A complete listing of those 


options is given in Appendix A. 


35 Input 

The input for ZTRAN is organized in a manner similar to that for ETRAN, 
and consists of three parts: (a) run parameters, (b) photon cross-section 
data, and (c) electron cross-section data. Items (a) and (b) are read in on 
logical unit 5 in card-image format. Item (c) is read in on logical unit 9, 


and consists of a data file generated by the pre-processing program DATAPAC. 


DATAPAC must be run with appropriate input parameters so as to include the 
materials and energies of interest. These input parameters are not discussed 
here, but are described in Ref. 2. 

A listing of the input from categories (a) and (b) is given in 
Appendix A. The run parameters include those for selecting the Monte Carlo 
model options and the parameters which determine the selection and size of the 
target layers, the type of incident radiation (electron or photons), the 
nature of initial spectrum (monoenergetic, continuous or discrete spectra), 
the type of initial angular distribution (mono-directional, isotropic, cosine- 
law), the position of source plane (external or internal), and the type of 


information and scoring regions desired. 


4, Output 

The print-out from ZTRAN includes an echo of much of the input data, some 
auxiliary intermediate information, and the various results including: 
(a) number and energy transmission and reflection coefficients for electrons 
and for photons crossing the target slab boundaries; (b) energy and charge 
deposition as a function of depth in the target; (c) spectra of electron and 
of photon fluence, as a function of depth, in the target; (d) spectra of 
energy absorbed in the target slabs; and (e) energy spectra, angular distribu- 
tions, and joint energy and angular distributions of electrons and of photons 
crossing slab boundaries in both forward and backward directions. 

The print-out includes sufficient headings so that it is largely self- 
explanatory. In lieu of a lengthy description of the ZTRAN output, we give in 
Appendix B the printout for a sample run, which was in fact used to solve one 


of the illustrative problems described in section 5. 


SIE Illustrative Applications of ZTRAN 

As a first example, we illustrate the application of ZTRAN to the 
calculation of depth-dose distributions from 60-MeV electron beams in 
composite media. Figure 1 compares the depth-dose curves in homogeneous 
water, homogeneous aluminum, and a composite target consisting of 2 g/cm? of 
aluminum followed by a semi-infinite water medium. Figures 2 and 3 give 
corresponding results for copper-water and lead-water combinations. 

A second example pertains to the use of dye-film dosimeters to determine 
the absorbed dose from low-energy electron beams. The assumed configuration, 
as indicated in figure 4, involves a 400-keV electron beam that passes through 
a l-mil titanium window and 10 cm of air, traverses a 2-mil nylon-based dye- 
film dosimeter, and is backscattered or stopped by a thick backing plate. The 
ZTRAN calculation produces, among other quantities, the energy deposited in 
the film dosimeter, and the charge deposited in the backing plate, per 
incident electron. Table 1 gives these results for various backing-plate 
materials. Results of this type could be used, for example, to calibrate the 


dye-film dosimeter or to calibrate the current from the accelerator. 
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Table l. 


Energy and charge deposition from a 400-keV electron beam incident 
on the composite target of figure 4. The results include: 


(a) The ratio D/J, of the absorbed-dose rate in the dye-film 
dosimeter to the electron current incident on the titanium 
window; 


(b) The "relative dose," defined as the dose rate D in the presence 
of the specified metal backing plate divided by the 
corresponding dose rate in the presence of a nylon backing 
plate; 


(c) The ratio J/Jyg of the current induced in the backing plate to 
the current incident on the titanium window; 


(d) The ratio D/J of the dose rate in the dye-film dosimeter to the 
current induced in the backing plate. 


Backing D/Jy Relative D/J 
Material (MeV g-1 cm2) Dose J/J (MeV g-1 cm2) 

nylon 4.56 1.00 0.847 5.38 

Al 5.39 1.18 0.765 7.05 

Fe 6.41 1.41 0.657 9.76 

Cu 6.61 1.45 0.638 10.4 

Ag Ta? 1.61 0.559 geil 

Au Sodls 1.79 0.464 17.6 

U 8.50 1.86 0.426 20.0 
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Fig. 4 


Comparison of depth-dose distributions from a 60-MeV electron beam. 
The dashed curves represent the distributions in homogeneous semi- 
infinite aluminum or water media. The histogram represents the 
distribution in a multi-layer medium consisting of 2 g/cm? of aluminum 
followed by a semi-infinite water layer. 


Same as figure 1, but for a copper-water combination. 
Same as figure 1, but for a lead-water combination. 


Target arrangement for the film dosimeter study. 
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APPENDIX A 


FORMAT FOR CARD IMAGE INPUT (LOGICAL UNIT 5) FOR ZTRAN 4C 


CROSS SECTION INPUT 


MSET 16 
LOOP TO ## MSET TIMES 
NSET,NTAB,NEDG 316 
MPAIR, C(MTAX(I),1I=1,NTAB) 1216 
CCEBCM,I),M=1,MTAXCI)),1I=1,NTAB) 6F12.6 (NEW CARD FOR EACH I) 
CCATBCM,1I),M=1,MTAXC1I)),1=1,NTAB) 6F12.6 (NEW CARD FOR EACH I) 
CCSURVB(M,1I),M=1,MTAXCI)),1I=1,NTAB) 6F12.6 (NEW CARD FOR EACH I) 
CPAIRBCM) ,M=1,MPAIR) 6F12.6 
LOOP TO # NEDG TIMES 
EDGK, PKEG, WK 3F12.6 
NUMK, NUMA 216 
CEKON) ,N=1,NUMK) 6F12.6 
CRKCN) ,N=1,NUMK) 6F12.6 
CEAUGCN),N=1,NUMA) 6F12.6 
CRAUGCN),N=1,NUMA) # 6F12.6 
DETOUR, EXCI1,AVGKE ## 3F12.6 
RUN PARAMETER INPUT 
IRNMAX 16 


THE PARAMETER IRNMAX GIVES THE NUMBER OF RUNS TO BE MADE WITH 
THE CROSS SECTIONS FOR THE GIVEN SET OF MATERIALS. FOR EACH 
OF THE RUNS, INPUT PARAMETERS MUST BE SUPPLIED AS GIVEN 
BELOW. INPUT DATA INDICATED BY (%) ARE TO BE INCLUDED ONLY 
UNDER CERTAIN OPTIONS CONTROLED BY THE VALUES OF OTHER INPUT 
PARAMETERS (SEE MEANING OF INPUT VARIABLES). 


COMMENT 72H 
NRUN, INC, ISTRG, INOK, IPHOT,IBAD, 
IDES,ITER,IPSC,IPUN,ITAPE 1116 


LMAX, ITMK, JMAX, IPMK, JPMAX, ISMK, 
JSMAX, TAMK, KMAX, IBMK, KPMAX, 
LMIX 1216 
IMAX,_LDEP,LZMAX,LEFLUX, LEFMAX, 
NCYC,LPFLUX,LPFMAX,NPCYC, 


NFCOS,LBIN 1116 
INRAN T12 
TIN, TCUT, TPCUT, TSAVE, ANGIN, ZIN GFA 25 
(%) JSPEC, ISAM, COMMENT 216,66H 
CX) CSPECINCJ),J=1,JSPEC) 6F12.5 
(%) CESPCJ),J=1,JSPEC) 6F12.5 
BNUM, XNUM, DLIM SRL225 
(BCL),L=1,LMAX) 6F12.5 
CMATNOCL),L=1,LMAX) 6112 
CLESPEC(L),L=1,LMAX) 1216 
CLPSPEC(L),L=1,LMAX) 1216 
C%) CTMARKCJ),J=1,JMAX) 6F112.5 
(%) CPMARKCJ),J=1,JPMAX) 6F12.5 
(*) CSMARKCJ),J=1,JSMAX) Ghlceo 
(%) CAMARKCK),K=1,KMAX) 6F12.5 
(%) CBMARKCK),K=1,KPMAX) 62.5 
(*) CBDCL),L=1,-LZMAX) 625 
(X) (CBFCL),L=1,-LEFMAX) GiRa2. 5 
(%) CFMARKCJ),J=1,JFMAX) Gir hee5 
(%) CBPFCL),L=1,-LPFMAX) GRL2 5 
(%) CFPMARK(CJ),J=1,JPFMAX) G2 /.5 
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MEANING OF INFUT VARIABLES FUR FRUGORAM CIRAN GC 


MSET 
NSET 
NTAB 
NEDG 
MPAIR 
MTAXCT) 
EB 


ATB 
SURVB 


PAIRB 


DETOUR 
EXCI1 


AVGKE 


IRNMAX 
COMMENT 
NRUN 
INC 


NUMBER OF MATERIALS FOR WHICH CROSS SECTIONS ARE TO BE 
READ. 

NUMBER OF DATAPAC SET FOR THE DESIRED MATERIAL THAT IS 
PRESENT ON THE INPUT FILE CUNIT 9). 

NUMBER OF ENERGY REGIONS FOR WHICH GAMMA RAY CROSS SECTION 
TABLES ARE GIVEN CEQUALS NUMBER OF ABSORPTION EDGES + 1). 
NUMBER OF SHELLS FOR WHICH CHARACTERISTIC X-RAY AND AUGER 
EMISSION WILL BE CONSIDERED. 

NUMBER OF ELEMENTS IN PAIRB ARRAY. 

NUMBER OF ELEMENTS IN I'TH TABLE OF EB ARRAY. 

LIST OF ENERGIES (MEV, DESCENDING ORDER) AT WHICH GAMMA 
RAY CROSS SECTIONS ARE TABULATED. 

TOTAL ATTENUATION COEFFICIENT (CMXX2/G). 

SUM OF COMPTON SCATTERING AND PAIR PRODUCTION COEFFICIENT 
DIVIDED BY TOTAL ATTENUATION COEFFICIENT CEQUALS 
PROBABILITY THAT AN INTERACTION WILL NOT RESULT IN 
PHOTOELECTRIC ABSORPTION). 

RATIO OF COMPTON SCATTERING COEFFICIENT TO SUM OF COMPTON 
SCATTERING AND PAIR PRODUCTION COEFFICIENTS CEQUALS 
PROBABILITY THAT A SCATTERING RATHER THAN A PAIR 
PRODUCTION EVENT OCCURS, EXCLUDING THE POSSIBILITY OF 
PHOTOELECTRIC ABSORPTION). 

SHELL BINDING ENERGY CIN MEV). 

SHELL PHOFOEFFECT EFFICIENCY. CEQUALS RATIO OF PHOTO- 
ELECTRIC CROSS SECTION FOR SHELL TO TOTAL PHOTOELECTRIC 
CROSS SECTION. THIS RATIO IS ASSUMED CONSTANT AND RELATES 
TO ENERGIES ABOVE THE K-SHELL BINDING ENERGY, I.E., WHERE 
ALL SHELLS CONTRIBUTE. AT LOWER ENERGIES, NON-CONTRIBUTING 
SHELLS ARE AUTOMATICALLY ACCOUNTED FOR.) 

SHELL FLUORESCENCE EFFICIENCY. 

ENERGIES CIN MEV) OF CHARACTERISTIC X-RAYS. 

CUMULATIVE PROBABILITY DISTRIBUTION FOR X-RAY ENERGIES. 
ENERGIES CIN MEV) OF AUGER ELECTRONS. 

Seer eae PROBABILITY DISTRIBUTION FOR AUGER ELECTRON 


THE PRESENT SCHEME OF X-RAY AND AUGER EMISSION IS NEITHER 
COMPLETE NOR GENERAL. FOR EXAMPLE, THE REARRANGEMENT OF 
THE ATOM FOLLOWING IONIZATION EVENTS IS NOT INCLUDED CE.G. 
L-SHELL EMISSION FOLLOWING K-SHELL FILLING, ETC.). IFA 
MATERIAL IS A COMPOUND OR MIXTURE, SHELLS BEYOND THE K- 
SHELL SHOULD NOT BE CONSIDERED. IN THIS CASE, THE K-ION- 
IZATION EVENTS CAN BE CALCULATED ONLY FOR THE ELEMENT WITH 
LARGEST Z. THE VARIABLES EDGK, WK, EK, RK, EAUG AND RAUG 
SHOULD BE GIVEN THE VALUES FOR THIS ELEMENT, AND PKEG 
GIVEN THE VALUE OF THE K-SHELL PHOTOEFFECT EFFICIENCY FOR 
THIS ELEMENT MULTIPLIED BY PROBABILITY THAT THE PHOTO- 
ELECTRIC ABSORPTION TAKES PLACE IN THE LARGEST-Z ELEMENT. 
THIS PROBABILITY IS 

(WXSIGPE)LARGEST-Z 7 SUMMATIONCWXSIGPE)D, 
WHERE W IS THE FRACTION BY WEIGHT AND SIGPE THE PHOTO- 
EFFECT CROSS SECTION IN BARNS/ATOM FOR AN ELEMENT. IN 
GENERAL, THIS PROBABILTY IS A FUNCTION OF PHOTON ENERGY. 
THEREFORE, AN AVERAGE OR EFFECTIVE PROBABILITY, SUITABLE 
TO THE PROBLEM CONSIDERED, SHOULD BE DETERMINED. 


RATIO OF PRACTICAL TO MEAN RESIDUAL ELECTRON RANGE. 

MEAN EXCITATION ENERGY CKEV) OF THE MEDIUM FOR COLLISION 
ENERGY-LOSS STRAGGLING. IF OMITTED, VALUE IS ESTIMATED 
INTERNALLY. 

MEAN KINETIC ENERGY CKEV) OF ATOMIC ELECTRONS IN THE 
MEDIUM, USED IN EVLUATION OF COLLISION ENERGY-LOSS 
STRAGGLING. IF OMITTED, VALUE IS ESTIMATED INTERNALLY. 


NUMBER OF RUNS TO BE MADE. 

ANY MESSAGE CONSISTING OF UP TO 72 HOLLERITH CHARACTERS. 
RUN NUMBER CUSED ONLY FOR IDENTIFICATION). 

1 CASCADE INITIATED BY ELECTRONS INCIDENT ON TARGET. 
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ISTRG 


INOK 


IPHOT 


IBAD 


IDES 


ITER 


IPSC 


IPUN 


ITAPE 


LMAX 
ITMK 


JMAX 
IPMK 


JPMAX 
ISMK 


JSMAX 
IAMK 


KMAX 


IBMK 


CASCADE INITIATED BY PHOTONS INCIDENT ON TARGET. 
ELECTRON ENERGY LOSS STRAGGLING TAKEN INTO ACCOUNT. 
ELECTRON ENERGY LOSS TREATED IN CONTINUOUS-SLOWING- 
DOWN APPROXIMATION. 

KNOCK-ON DELTA-RAY HISTORIES ARE FOLLOWED. 

KNOCK-ON DELTA-RAY HISTORIES ARE NOT FOLLOWED. 

HISTORIES OF SECONDARY PHOTONS (CBREMSSTRAHLUNG, 

CHARACTERISTIC X-RAYS, ANNIHILATION RADIATION) ARE 

FOLLOWED. 

HISTORIES OF SECONDARY PHOTONS ARE NOT FOLLOWED. 

INTRINSIC EMISSION ANGLE FOR BREMSSTRAHLUNG SAMPLED 

FROM DETAILED CROSS SECTION TABLES. 

INTRINSIC EMISSION ANGLE SAMPLED WITH USE OF SIMPLE 

APPROXIMATE FORMULA. 

ore TORS OF PHOTON-PRODUCED SECONDARY ELECTRONS ARE 

2 THESE HISTORIES ARE NOT FOLLOWED. 

pee TNES FATE OF ELECTRONS WITH ENERGIES BETWEEN TCUT 

I HISTORY IS CONTINUED ONLY IF ELECTRON CAN REACH A MAJOR 

BOUNDARY, B. 
2 HISTORY IS CONTINUED ONLY IF ELECTRON CAN REACH AN 
ENERGY-DEPOSITION SUB-BOUNDARY, BD. 

CONTROLS THE SCORING OF EMERGENT PHOTON DISTRIBUTIONS AND 

INTERNAL PHOTON FLUX. 

1 SCORE ALL PHOTONS. 

2 SCORE PRIMARY PHOTONS AND "PHIST DESCENDANTS" CANNIHIL- 
ATION QUANTA OR X RAYS RESULTING FROM PHOTON INTER- 
ACTIONS WITHOUT CONTROL PASSING FROM SUBROUTINE PHIST). 
SCORE BREMSSTRAHLUNG PHOTONS AND "PHIST DESCENDANTS". 
SCORE K X-RAYS AND "PHIST DESCENDANTS". 

SCORE ANNIHILATION QUANTA AND "PHIST DESCENDANTS". 

SOME RESULTS PUNCHED OUT ON CARDS. IT IS ASSUMED THAT 
ZTRAN IS SUPPLIED WITH A SUBROUTINE PUNOUT WHICH 
DETERMINES THE PUNCHING FORMAT AND THE TYPE OF 
INFORMATION INCLUDED. 

NO RESULTS PUNCHED ON CARDS. 

PROVISION IS MADE FOR WRITING RESULTS ON TAPE CUNIT 11) 
NO SPECIAL TAPE WRITTEN. 

UMBER OF SLABS MAKING UP COMPOSITE TARGET. 

ENERGY SPACING FOR TMARK GRID LINEAR FROM TIN TO ZERO. 
THE SMALLEST VALUE OF TMARK = TCUT. THE ENERGY 
DECREMENT = TIN/JMAX EXCEPT FOR THE LAST BIN. 
JMAX MAY BE REDUCED FROM ITS ORIGINAL INPUT VALUE 
DEPENDING ON THE VALUE OF TCUT. 
2 ENERGY SPACING FOR TMARK GRID LINEAR FROM TIN TO TCUT. 
3 ees FOR TMARK LOGARITHMICALLY SPACED FROM TIN 
0 TCUT. 
4 ENERGY SPACING FOR TMARK GRID READ IN AS TABLE. THE 
SMALLEST VALUE OF TMARK=TCUT IS INSERTED, AND JMAX MAY 
BE REDUCED. 

NUMBER OF ENERGY BINS (DEFINED BY TMARK ARRAY) FOR 

CLASSIFYING EMERGENT ELECTRONS. 

SAME SET OF OPTIONS AS FOR ITMK, BUT DEFINED WITH 

REFERENCE TO PMARK,TPCUT AND JPMAX INSTEAD OF TMARK, 

TCUT AND JMAX. 

NUMBER OF ENERGY BINS (DEFINED BY PMARK ARRAY) FOR 

CLASSIFYING EMERGENT PHOTONS. 

SAME SET OF OPTIONS AS FOR ITMK, BUT DEFINED WITH 

REFERENCE TO SMARK, TCUT CINC=1) OR TPCUT CINC=2) AND 

JSMAX INSTEAD OF TMARK, TCUT AND JMAX. 

NUMBER OF ENERGY BINS (DEFINED BY SMARK ARRAY) FOR 

CLASSIFYING ABSORBED ENERGY. 

1 ANGULAR GRID FOR AMARK LINEAR FROM 0 TO 90 DEGREES 

CANGLE INCREMEMT = 90/KMAX DEGREES). 

2 ANGULAR GRID FOR AMARK READ IN AS TABLE. 

NUMBER OF ANGULAR BINS FOR CLASSIFYING THE TRANSMITTED 

AND REFLECTED ELECTRONS IN REGARD TO THEIR OBLIQUITY 

WITH RESPECT TO THE NORMAL TO THE TARGET. 

SAME OPTIONS AS FOR IAMK, BUT DEFINED WITH REFERENCE TO 


RPNe NFA 


mm NY FN 


Mn pwW 


Sone N 
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KPMAX 
LMIX 


IMAX 


LDEP 


LZMAX 


LEFLUX 
LEFMAX 


NCYC 


LPFLUX 
LPFMAX 


NPCYC 


NFCOS 


LBIN 
INRAN 


TIN 


TCUT 


TPCUT 


TSAVE 


BMARK AND KPMAX INSTEAD OF AMARK AND KMAX. 

NUMBER OF BINS FOR CLASSIFYING THE TRANSMITTED AND 
REFLECTED PHOTONS IN REGARD TO THEIR OBLIQUITY. 

INDEX OF LAST SLAB IN WHICH PHOTONS ARE FOLLOWED CUSED FOR 
GENERATING DESCENDANT ELECTRON SOURCES). IF LMIX = 0, 
LMIX IS SET EQUAL TO LMAX. 

NUMBER OF HISTORIES OF PRIMARY PARTICLES TO BE GENERATED. 
THESE PRIMARIES MAY BE EITHER ELECTRONS OR PHOTONS, 
DEPENDING ON THE VALUE OF THE PARAMETER INC. 

FOR LDEP GREATER THAN 0, ENERGY AND CHARGE DEPOSITION ARE 
CALCULATED FOR THE REGION OCCUPIED BY SLABS 1 THROUGH LDEP 
FOR LDEP LESS THAN 0, ENERGY AND CHARGE DEPOSITION ARE 

at amy FOR THE REGION OCCUPIED BY THE LDEP-TH SLAB 

FOR LDEP=0, ENERGY AND CHARGE DEPOSITION ARE NOT CALCULATE 
AND PROGRAM SETS ITER=1. 

FOR LZMAX GREATER THAN 0, THE NUMBER OF EQUAL DEPTH- 
REGIONS INTO WHICH THE ENERGY AND CHARGE DEPOSITION REGION 
DEFINED BY LDEP IS DIVIDED. FOR LZMAX LESS THAN 0, -LZMAX 
Beer ane CHARGE DEPOSITION DEPTH-REGION BOUNDARIES ARE 
seetane TO LDEP, BUT DEFINING THE ELECTRON FLUX SCORING 
FOR LEFMAX GREATER THAN 0, NUMBER OF EQUAL SUB-REGIONS 
INTO WHICH THE ELECTRON FLUX REGION, DEFINED BY LEFLUX, IS 
DIVIDED. FOR LEFMAX LESS THAN 0, ~LEFMAX ELECTRON FLUX 
SUB-REGION BOUNDARIES ARE READ IN. THE DIRECTIONAL FLUX 
DISTRIBUTIONS ARE RECORDED FOR EACH SUB-REGION. 

1 OR GREATER. THE BOTTOM ENERGIES OF THE ENERGY CLASSIFI- 


CATION BINS FOR THE ELECTRON FLUX ARE EQUAL TO 
TINX(€2*XXC(-N/NCYC)). TCUT IS ADDED AS THE BOTTOM ENERGY 
INSURE BESi: 

OR SMALLER. THE FIRST DIGITS OF THIS NUMBER 

SUPPLY THE DESIRED NUMBER, JFMAX, OF ENERGY BINS FOR 
THE ELECTRON FLUX CALCULATION. THE LAST DIGIT ALLOWS 
Fee ENERGY CLASSIFICATION OPTIONS SIMILAR TO THOSE FOR 


TMK. 

ae TO LDEP, BUT DEFINING THE PHOTON FLUX SCORING 
SIMILAR TO LEFMAX BUT DETERMINES LPFMAX PHOTON FLUX SUB- 
REGION BOUNDARIES. FLUX DISTRIBUTIONS OF POSITIVELY AND 
ME CROH on. DIRECTED PHOTONS ARE RECORDED FOR EACH SUB- 
SAME OPTIONS AS FOR NCYC, BUT DEFINES ENERGY 
CLASSIFICATION BINS FOR THE PHOTON FLUX, JPFMAX IN NUMBER. 
TPCUT IS ADDED AS THE BOTTOM ENERGY. 
NUMBER OF EQUAL SOLID ANGLES INTO WHICH FULL SPACE IS 
DIVIDED FOR THE CLASSIFICATION OF THE ELECTRON DIRECTIONAL 
FLUX DISTRIBUTIONS. 

INDEX OF SLAB IN WHICH INITIAL SOURCE IS DEFINED. IF 
LBIN = 0, LBIN IS SET EQUAL TO l. 

INITIAL RANDOM NUMBER CODD INTEGER).ON IBM 7094 AND 
UNIVAC 1108 AN OCTAL 12-DIGIT NUMBER WHOSE FIRST DIGIT IS 
NO eeeER THAN 3. FOR THE IBM 360 SERIES, A DECIMAL 

NTEGER. 

ENERGY CIN MEV) OF PRIMARY PARTICLES CELECTRONS OR 
PHOTONS) INCIDENT ON TARGET. WHEN TIN IS NEGATIVE, 
INCIDENT SPECTRUM IS NOT MONOENERGETIC, AND PROGRAM 
EXPECTS INPUT SPECTRUM CJSPEC, SPECIN, ESP). 
CUT-OFF ENERGY CIN MEV) AT WHICH ELECTRON HISTORIES 
ARE TERMINATED, A FINAL ADJUSTMENT IS MADE PERTAINING 
TO THE CALCULATION OF ENERGY AND CHARGE DEPOSITION. THIS 
ADJUSTMENT IS SIMILAR TO THAT DESCRIBED IN THE COMMENT 

ON THE PARAMETER TSAVE GIVEN BELOW. 
CUT-OFF ENERGY CIN MEV) AT WHICH PHOTON HISTORIES 
ARE TERMINATED. UPON TERMINATION THE RESIDUAL ENERGY OF 
THE PHOTON IS ASSUMED TO ESCAPE THE TARGET. 
A CUT-OFF ENERGY CIN MEV) GREATER THAN OR EQUAL TO TCUT, 
WHICH BECOMES OPERATIONAL WHEN THE ELECTRON IS TRAPPED, 
I.E. WHEN THE RESIDUAL RANGE IS SMALLER THAN THE DISTANCE 
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ANGIN 


ZIN 


JSPEC 


ISAM 


SPECIN 


ESP 
BNUM 


XNUM 


DLIM 


TO THE NEAREST BOUNDARY (SEE ITER). THE HISTORIES OF TRAP 
ELECTRONS WITH ENERGIES LESS THAN TSAVE ARE TERMINATED. 
IF SECONDARY ELECTRONS ARE PRODUCED WITH ENERGIES LESS 
THAN TSAVE AND IN A TRAPPED CONDITION, THEIR HISTORIES ARE 
NOT FOLLOWED. WHEN HISTORIES ARE EITHER TERMINATED OR NOT 
FOLLOWED, A FINAL ADJUSTMENT IS MADE FOR THE PURPOSE OF 
SCORING ENERGY AND CHARGE DEPOSITION. THE ELECTRON IS 
MADE TO TRAVERSE A FINAL SHORT STRAIGHT PATH IN THE 
DIRECTION REACHED UPON TERMINATION COR AT BIRTH IN THE 
CASE OF SECONDARY ELECTRONS). THE LENGTH OF THE FINAL 
PATH SEGMENT IS TAKEN TO BE EQUAL TO THE RESIDUAL 
PRACTICAL RANGE CRESIDUAL MEAN RANGE MULTIPLIED BY A 
DETOUR FACTOR WHICH IS SUPPLIED AS AN INPUT DATUM), AND 
THE RESIDUAL ENERGY AND THE CHARGE OF THE ELECTRON ARE 
ALLOWED TO BE DEPOSITED SOMEWHERE AT RANDOM ALONG THE 
FINAL PATH SEGMENT. 
FOR LEFLUX NOT 0 CELECTRON FLUX DISTRIBUTIONS CALCULATED), 
EE eR OeR An SETS TSAVE=TCUT, REGARDLESS OF THE INPUT VALUE 
90.0 OR SMALLER, FOR MONODIRECTIONAL SOURCE. ANGIN IS THE 
ANGLE CIN DEGREES) OF INCIDENCE WITH RESPECT TO THE 
NORMAL TO THE TARGET. 
100.0 + ANGLE CIN DEGREES), FOR AN ISOTROPIC SOURCE 
BETWEEN O AND ANGLE. 
200.0 + ANGLE CIN DEGREES), FOR A COSINE LAW SOURCE 
BETWEEN 0 AND ANGLE. 
SOURCE PLANE POSITION, IN G/CM2. IF ZIN = 0, ZIN IS SET 
EQUAL TO BCLBIN+1). IF ZIN = LBIN = 0, ZIN IS SET EQUAL 
TO O AND LBIN IS SET EQUAL TO l. : 
NUMBER OF END POINTS OF THE ENERGY BINS IN THE SPECTRUM OF 
THE INCIDENT RADIATION. JSPEC, ISAM, SPECIN AND ESP ARE 
READ IN ONLY WHEN TIN IS NEGATIVE. 
1 IN SAMPLING FROM THE SOURCE SPECTRUM, INTERPOLATION 
LINEAR IN ENERGY IS USED. 
2 IN SAMPLING FROM THE SOURCE SPECTRUM, INTERPOLATION 
LOGARITHMIC IN ENERGY IS USED. 
3 SRE ERC C TRUM IS SAMPLED FROM A SET OF DISCRETE 
ENERGIES. 
FOR ISAM=1 AND 2, CUMULATIVE PROBABILITY DISTRIBUTION FOR 
THE SPECTRUM OF INCIDENT RADIATION. SPECINC1)=1.0 AND 
SPECINCJSPEC)=0.0. FOR ISAM=3, RELATIVE PROBABILITIES 
OF THE DISCRETE LINES COMPRISING THE INCIDENT SPECTRUM. 
THESE PROBABILITIES ARE ALTERED SLIGHTLY BY THE PROGRAM 
SO THAT INTEGRAL NUMBERS OF HISTORIES ARE CALCULATED FOR 
EACH LINE. THE NUMBERS OF HISTORIES AND THE ALTERED 
RELATIVE PROBABILITIES ARE PRINTED OUT FOR REFERENCE. 
ENERGIES CIN MEV) IN DESCENDING ORDER ASSOCIATED WITH THE 
DISTRIBUTION SPECIN. 
NOMINAL NUMBER OF BREMSSTRAHLUNG PHOTONS (WITH ENERGIES 
GREATER THAN TPCUT) THAT ONE WANTS TO SAMPLE PER ELECTRON 
HISTORY. USING THE BREMSSTRAHLUNG CROSS SECTION AND 
ASSUMING THE CONTINUOUS SLOWING-DOWN APPROXIMATION FOR THE 
ELECTRON ENERGY LOSS BETWEEN TIN AND TCUT, THE PROGRAM 
CALCULATES A FACTOR CALLED SCALE. THE BREMSSTRAHLUNG 
CROSS SECTIONS ARE THEN MULTIPLIED BY SCALE, AND THE 
SAMPLED PHOTONS ARE GIVEN THE STATISTICAL WEIGHT 
1.0/7SCALE. THIS ADJUSTMENT IS MADE ON THE BASIS OF 
PRODUCTION PROBABILITIES PERTAINING ONLY TO THE MATERIAL 
OF THE FIRST SLAB AND DOES NOT TAKE INTO ACCOUNT PASSAGE 
OF THE ELECTRONS INTO SUBSEQUENT SLABS. DUE TO THIS 
EFFECT, THE OCCURRENCE OF ENERGY-LOSS STRAGGLING, 
AND THE ESCAPE OF ELECTRONS FROM THE TARGET, 
THE ACTUAL NUMBER OF PHOTONS SAMPLED IN A RUN WILL DIFFER 
SOMEWHAT FROM THE EXPECTED NOMINAL NUMBER. THE PROGRAM 
DOES NOT ALLOW SCALE TO HAVE A VALUE SMALLER THAN 1.0. 
THE VALUE SCALE = 1.0 CAN BE OBTAINED BY SETTING BNUM=0.0. 
PARAMETER ANALOGOUS TO BNUM, BUT FOR K-SHELL IONIZATION BY 
ELECTRON IMPACT. THE NOMINAL NUMBER OF EMITTED K X-RAYS 
WILL BE WKXXNUM. 
PROBABILITY OF FOLLOWING THE HISTORY OF A SECONDARY 
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B 
MATNO 


LESPECCL) 


LPSPECCL) 


TMARK 


PMARK 


SMARK 


AMARK 


BMARK 


BD 


BF 


FMARK 


BPF 


FPMARK 


ELECTRON PRODUCED BY A PHOTON. IF THE ENERGY OF THE 

NEWLY BORN ELECTRON IS GREATER THAN TSAVE, OR IF IT IS 

BETWEEN TCUT AND TSAVE AND THE RESIDUAL RANGE IS LONG 

ENOUGH FOR THE NEAREST BOUNDARY TO BE REACHED CSEE ITER), 

THE HISTORY IS FOLLOWED WITH PROBABILITY DLIM. THE 

STATISTICAL WEIGHT ASSOCIATED WITH THE HISTORY IS 

MULTIPLIED BY 1.0/DLIM. REGARDLESS OF THE INPUT VALUE OF 

DLIM, THE PROGRAM ADJUSTS ITS VALUE, IF NECESSARY, SO THAT 

DLIM IS GREATER THAN 1.0/SCALE. THE VALUE DLIM=1.0/SCALE 

CAN BE OBTAINED BY SETTING DLIM=0.0. 

THICKNESSES OF SLABS DEFINING TARGET, IN G/CMxXX2 AND IN 

ASCENDING ORDER. 

INDICES SPECIFYING WHICH OF THE MSET CROSS-SECTION DATA 

SETS PERTAIN TO THE SLABS B. THE NUMBERING CORRESPONDS TO 

THE ORDER IN WHICH THESE DATA SETS WERE READ. 

PERTAINS TO JOINT ANGULAR AND ENERGY DISTRIBUTION OF 

ELECTRONS EMERGING FROM L*TH TARGET SLAB. 

0 NO DISTRIBUTION CALCULATED. 

1 DISTRIBUTION CALCULATED FOR TRANSMITTED ELECTRONS. 

2 DISTRIBUTION CALCULATED FOR REFLECTED ELECTRONS. 

3 DISTRIBUTION CALCULATED FOR TRANSMITTED AND REFLECTED 

ELECTRONS. 

SIMILAR TO LESPEC(L), BUT PERTAINS TO JOINT ANGULAR AND 

eee | DISTRIBUTION OF PHOTONS EMERGING FROM L'TH TARGET 

SLAB. 

BOTTOM ENERGY (MEV, DESCENDING ORDER) OF ENERGY BINS 

FOR CLASSIFYING SPECTRA OF TRANSMITTED AND REFLECTED 

ELECTRONS. READ IN AS TABLE ONLY WHEN ITMK=4. 

BOTTOM ENERGY (MEV, DESCENDING ORDER) OF ENERGY BINS 

FOR CLASSIFYING SPECTRA OF TRANSMITTED AND REFLECTED 

PHOTONS. READ IN AS TABLE ONLY WHEN IPMK=4. 

BOTTOM ENERGY (MEV, DESCENDING ORDER) OF ENERGY BINS 

THAT ARE THE COMPLIMENT OF THOSE USED TO CLASSIFY ABSORBED 

ENERGY (I.E., TIN OR ESPC1) MINUS SMARKCJ)). READ IN AS 

TABLE ONLY WHEN ISMK=4. 

TOP OF ELECTRON OBLIQUITY BINS CIN DEGREES), IN ASCENDING 

aha TO 90 DEGREES. TO BE READ IN AS TABLE ONLY WHEN 

AMK=2. 

TOP OF PHOTON OBLIQUITY BINS CIN DEGREES), IN ASCENDING 

ae 90 DEGREES. TO BE READ IN AS TABLE ONLY WHEN 
MK=2. 

ENERGY AND CHARGE DEPOSITION SUB-REGION BOUNDARIES CIN 

ASCENDING ORDER) IN G/CMxXxX2 CFOR IUNT=2). TO BE READ IN 

AS TABLE ONLY WHEN LZMAX IS LESS THAN 0O. 

ELECTRON FLUX SUB-REGION BOUNDARIES CIN ASCENDING ORDER) 

IN G/CMXX¥2 (FOR IUNT=2). TO BE READ IN AS TABLE ONLY WHEN 

LEFMAX IS LESS THAN 0. 

BOTTOM ENERGY (MEV, DESCENDING ORDER) OF ENERGY BINS FOR 

CLASSIFYING ELECTRON FLUX SPECTRA. READ IN AS TABLE ONLY 

FOR NCYC LESS THAN 0 WITH LAST DIGIT EQUAL TO 4G. 

PHOTON FLUX SUB-REGION BOUNDARIES CIN ASCENDING ORDER) IN 

G/CM¥X2 (FOR IUNT=2). TO BE READ IN ONLY WHEN LPFMAX IS 

LESS THAN 0. 

BOTTOM ENERGY (MEV, DESCENDING ORDER) OF ENERGY BINS FOR 

CLASSIFYING PHOTON FLUX SPECTRA. READ IN AS TABLE ONLY 

FOR NPCYC LESS THAN 0 WITH LAST DIGIT EQUAL TO 4. 
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BRIEF DESCRIPTION OF SUBROUTINES USED IN PROGRAM ZTRAN 4C 


SUBROUTINE RANDA 
GENERATES PSEUDO-RANDOM NUMBERS. 


SUBROUTINE QPOL 
PERFORMS PARABOLIC INTERPOLATION. 


SUBROUTINE CLASS 
FROM ELECTRON ENERGY, CALCULATES INDEX REQUIRED FOR TABLE LOOK-UP OR 
EVALUATION OF QUANTITIES DEPENDING ON ELECTRON ENERGY. 


SUBROUTINE CROSS 


CALCULATES ENERGY INDEX REQUIRED FOR TABLE LOOK-UP OF PHOTON CROSS 
SECTIONS. 


SUBROUTINE XINPUT 
READS INTO MEMORY VARIOUS CROSS SECTIONS AND MULTIPLE SCATTERING 
DISTRIBUTIONS. INPUT COMES FROM CARDS AND FROM DATAPAC TAPE. 


SUBROUTINE XPREP 
PROCESSES INPUT DATA TO FACILITATE QUICK AND EASY TABLE LOOK-UP. 


SUBROUTINE PTAB 
PREPARES TABLES OF PHOTON CROSS SECTIONS NEEDED FOR TABLE LOOK-UPS. 


SUBROUTINE INPUT 
READS INTO MEMORY THE VARIOUS RUN PARAMETERS. 


SUBROUTINE START 
COMPUTES INITIAL ELECTRON RANGE. CHECKS CONSISTENCY OF VARIOUS RUN 
Repo pes AND HALTS RUNNING OF PROGRAM IF INCONSISTENCIES ARE 


SUBROUTINE FLIST 
GENERATES, OR READS IN AS INPUT, ENERGY LISTS REQUIRED FOR 
CLASSIFYING ELECTRON AND PHOTON FLUX DATA. 


SUBROUTINE ELIST 
GENERATES, OR READS IN AS INPUT, ENERGY LISTS REQUIRED FOR 
CLASSIFYING ELECTRON AND PHOTON OUTPUT DATA. 


SUBROUTINE ALIST 
GENERATES, OR READS AS INPUT, ANGLE LISTS REQUIRED FOR 
CLASSIFYING ELECTRON AND PHOTON OUTPUT DATA. 


SUBROUTINE PREP 
DERIVES FROM RUN PARAMETERS AND CROSS SECTION DATA VARIOUS ADDITIONAL 
AUXILIARY TABLES AND PARAMETERS. ALSO CLEARS OUTPUT STORAGE ARRAYS. 


SUBROUTINE SECP 
CALCULATES AUXILIARY TABLES PERTAINING TO THE PRODUCTION OF ENERGETIC 
KNOGK-ON ELECTRONS. 


SUBROUTNE CLASSP 
FROM PHOTON ENERGY, CALCULATES INDEX REQUIRED FOR TABLE LOOK-UP OR 
EVALUATION OF QUANTITIES DEPENDING ON ELECTRON ENERGY. 


SUBROUTINE CRISS 
COMPUTES PARAMETERS NEEDED FOR SPEEDING UP TABLE LOOK-UP IN 
SUBROUTINE CROSS 


SUBROUTINE HIST 

SETS UP INITIAL CONDITIONS FOR PRIMARY AND SECONDARY ELECTRONS IN 
CASE OF AN ELECTRON-INITIATED CASCADE, OR FOR PRIMARY PHOTONS AND 
SECONDARY ELECTRONS IN CASE OF A PHOTON-INITIATED CASCADE. CALLS 
THE VARIOUS ROUTINES REQUIRED FOR THE GENERATION OF PRIMARY AND 
SECONDARY ELECTRON HISTORIES. KEEPS TRACK OF THE TOTAL ENERGY 
DEPOSITED IN THE TARGET BY THE PARTICLE INITIATING THE CASCADE AND 
ALL ITS DESCENDANTS. 
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SUBROUTINE EHIST 

GENERATES A SINGLE ELECTRON HISTORY CPRIMARY OR SECONDARY) AND 

CALLS THE SUBROUTINES REQUIRED TO CALCULATE THE DESIRED OUTPUT INFOR- 
MATION CTRANSMISSION, REFLECTION, ENERGY AND CHARGE DEPOSITON, FLUX). 


SUBROUTINE LOSS 
SAMPLES ELECTRON ENERGY LOSS FROM LANDAU DISTRIBUTION CWITH BLUNCK- 
LEISEGANG CORRECTION). 


SUBROUTINE MULT 
SAMPLES ANGULAR DEFLECTION DUE TO MULTIPLE ELASTIC SCATTERING FROM 
GOUDSMIT-SAUNDERSON DISTRIBUTION. 


SUBROUTINE KNOCK 

SAMPLES INITIAL ENERGY AND DIRECTION OF KNOCK-ON ELECTRON ACCORDING 
TO MOLLER CROSS SECTION. IF THE HISTORY OF THE SECONDARY PARTICLE IS 
NOT TO BE FOLLOWED, DISPOSES OF THE RESIDUAL ENERGY FOR THE PURPOSE 
OF SCORING ENERGY DEPOSITION. 


SUBROUTINE BREMS 

CASE A CBNUM=0.0, BREMSSTRAHLUNG PHOTONS SAMPLED IN 

ACCORDANCE WITH NATURAL OCCURRENCE). SAMPLES NUMBER OF PHOTONS 
PRODUCED IN STEP FROM POISSON DISTRIBUTION WITH MEAN GIVEN BY BREMS- 
STRAHLUNG CROSS SECTION. SCORES INTERNAL BREMSSTRAHLUNG EFFICIENCY, 
SUPPLIES BREMSSTRAHLUNG ENERGY LOSS FOR SUBROUTINE EHIST, SAMPLES 
INITIAL DIRECTIONS, ENERGIES AND POSITIONS FOR BREMSSTRAHLUNG 
PHOTONS. IF THE SECONDARY PHOTON HISTORIES ARE TO BE FOLLOWED, 
CALLS REQUIRED SUBROUTINES TO ACCOMPLISH THIS. 

CASE B (CBNUM GREATER THAN 0, ARTIFICIALLY INCREASED NUMBER OF PHOTONS 
SAMPLED). AS IN CASE A, SUPPLIES BREMSSTRAHLUNG ENERGY LOSS FOR 
SUBROUTINE EHIST. GENERATES INITIAL CONDITIONS FOR SCALE TIMES THE 
MEAN NUMBER OF BREMSSTRAHLUNG PHOTONS WITH ENERGIES ABOVE TPCUT. 
THESE PHOTONS ARE SUPPLIED WITH APPROPRIATE STATISTICAL WEIGHTS AND 
FOLLOWED IN ORDER TO OBTAIN INFORMATION ABOUT THE INTERNAL BREMS~ 
STRAHLUNG EFFICIENCY, THICK-TARGET BREMSSTRAHLUNG AND THE SECONDARY 
ELECTRONS PRODUCED BY THE BREMSSTRAHLUNG. 


SUBROUTINE PHOTS 
SAMPLES ENERGY OF PHOTON GENERATED IN BREMSSTRAHLUNG EVENT. 


SUBROUTINE PINT 
SCORES INTERNAL BREMSSTRAHLUNG EFFICIENCY AND CALLS SUBROUTINE 
PHIST FOR THE GENERATION OF PHOTON TRAJECTORIES. 


SUBROUTINE KXRAY 

SIMILAR IN ITS TASKS TQ SUBROUTINE BREMS, BUT PERTAINS TO EMISSION OF 
K X-RAYS DUE TO ELECTRON-IMPACT IONIZATION RATHER THAN TO BREMS- 
STRAHLUNG PHOTONS FROM THE CONTINUOUS SPECTRUM. 


SUBROUTINE PSTART 

SAMPLES POSITION WHERE SECONDARY PHOTONS ARE BORN AND ESTIMATES 
ELECTRON DIRECTION AT THAT POINT. REJECTS PHOTONS BORN OUTSIDE 
THE MEDIUM. 


SUBROUTINE ESCORE 
SCORES DATA PERTINENT TO THE TRANSMISSION AND REFLECTION OF 
ELECTRONS. 


SUBROUTINE ANGLE 

SAMPLES ANGULAR MULTIPLE SCATTERING DEFLECTION FROM EXPONENTIAL 
APPROXIMATION TO GOUDSMIT-SAUNDERSON DISTRIBUTION, TO BE USED IN 
PARTIAL SUBSTEPS TO A BOUNDARY. 


SUBROUTINE SPILL 
SCORES ENERGY DEPOSITION BY ELECTRONS IN THE DEFINED ENERGY 
DEPOSITION SUB-REGIONS. 


SUBROUTINE RAID 
SIMILAR TO SPILL, BUT SCORES CHARGE DEPOSITION. 
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SUBROUTINE FORN 


SCORES DIRECTIONAL FLUX SPECTRUM OF ELECTRONS IN THE DEFINED ELECTRON 
FLUX SUB-REGIONS. 


SUBROUTINE PHIST 

MAIN ROUTINE FOR SCORING PHOTON HISTORIES. TAKES INTO ACCOUNT PHOTO- 
ELECTRIC ABSORPTION, COMPTON SCATTERING AND PAIR PRODUCTION, AS WELL 
AS THE SECONDARY ELECTRONS AND PHOTONS RESULTING FROM THESE EVENTS. 
THE INITIAL CONDITIONS FOR SECONDARY ELECTRONS ARE EXAMINED BY 
SUBROUTINE SAVE TO DETERMINE WHETHER THEIR HISTORIES SHOULD BE 
FOLLOWED AT ALL. IF THE ANSWER IS YES, THE INITIAL CONDITIONS ARE 
STORED FOR LATER USE BY SUBROUTINE HIST, BUT ONLY IN A FRACTION DLIM 
OF ALL CASES. IF THE ANSWER IS NO, THE CHARGE AND RESIDUAL PARTICLE 
ENERGY ARE TAKEN INTO ACCOUNT FOR THE PURPOSE OF SCORING ENERGY 

AND CHARGE DEPOSITION. THE POSSIBILITY OF THE EMISSION OF PAIR 
ANNIHILATION AT THE END OF A POSITRON TRACK IS ALSO TAKEN INTO 
ACCOUNT, AND THE RESULTING PHOTONS ARE FOLLOWED IN SUBROUTINE PHIST. 
IN THE PRESENT VERSION, THE POSSIBILITY OF ANNIHILATION 

IN FLIGHT IS DISREGARDED. IF FLUORESCENCE PHOTONS ARE PRODUCED 
AFTER A PHOTOELECTRIC ABSORPTION EVENT, THEIR HISTORIES ARE 

ALSO FOLLOWED IN SUBROUTINE PHIST. 


SUBROUTINE PATH 
DETERMINES PATHLENGTH THROUGH MULTI-LAYERED TARGET TO NEXT PHOTON 
INTERACTION. SIGNALS END OF PHOTON HISTORY DUE TO ESCAPE. 


SUBROUTINE PSCORE 

SCORES PHOTON TRANSMISSION AND REFLECTION WITH THE USE OF AN WEIGHT 
FACTOR EQUAL TO THE PROBABILITY THAT THE PHOTON, SUBSEQUENT TO THE 
COLLISION UNDER CONSIDERATION, WILL REACH THE BOUNDARY OF INTEREST 
WITHOUT FURTHER INTERACTION IN THE MEDIUM. FOR THE SCORING OF 
ABSORBED SPECTRA, HOWEVER, THE PHOTON HISTORY MUST CROSS A BOUNDARY 
TO COUNT AS EMERGENT ENERGY. 


SUBROUTINE PORN 
SCORES FLUX SPECTRUM OF POSITIVELY AND NEGATIVELY DIRECTED PHOTONS 
IN THE DEFINED PHOTON FLUX SUB-REGIONS. 


SUBROUTINE BANGLE 

SAMPLES THE INTRINSIC BREMSSTRAHLUNG PHOTON EMISSION ANGLE CWITH 
RESPECT TO THE DIRECTION OF THE ELECTRON) FROM TABULATED DISTRIBUTION 
(DERIVED FROM BETHE-HEITLER THEORY CROSS SECTIONS). 


SUBROUTINE LANGLE 

SIMILAR TO BANGLE, BUT USES SIMPLE APPROXIMATION TO BETHE-HEITLER 
CROSS SECTION. IS ALSO USED TO SAMPLE THE INTRINSIC EMISSION 
ANGLE FOR PAIR ELECTONS CWITH RESPECT TO THE PHOTON DIRECTION). 


SUBROUTINE RANGLE 
COMPUTES RANDOM DIRECTION COSINE DISTRIBUTED UNIFORMLY BETWEEN -1 
AND +1. 


SUBROUTINE PEAL 

SAMPLES DIRECTION OF PHOTO-ELECTRON FROM TABULATED DISTRIBUTIONS 
(BASED ON SAUTER CROSS SECTION AT HIGH ENERGIES AND ON FISCHER CROSS 
SECTION AT LOW ENERGIES). 


SUBROUTINE SPLIT 
SAMPLES ENERGY OF PAIR ELECTRON FROM TABULATED DISTRIBUTIONS (BASED 
ON HOUGH CROSS SECTION). 


SUBROUTINE SCATT 

SAMPLES ENERGY OF COMPTON RECOIL ELECTRON FROM KLEIN-NISHINA 
DISTRIBUTION, AND COMPUTES CORRESPONDING ANGULAR DEFLECTION OF THE 
PHOTON. 


SUBROUTINE FOLD 


COMPOUNDS INITIAL OBLIQUITY WITH SAMPLED DEFLECTION TO GIVE FINAL 
OBLIQUITY IN FIXED COORDINATE SYSTEM. 
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SUBROUTINE SAVE 4 

MAKES THE DECISION WHETHER A HISTORY SHOULD BE CONTINUED, AND 

WHETHER THE INITIAL CONDITIONS FOR A NEW HISTORY SHOULD BE SAVED FOR 
LATER USE. HISTORIES ARE CONTINUED OR SAVED WHEN THE ELECTRON 

ENERGY IS GREATER THAN TSAVE, OR WHEN THE ENERGY IS SMALLER THAN 
TSAVE BUT GREATER THAN TCUT AND THE RESIDUAL RANGE IS LONG ENOUGH FOR 
THE ELECTRON TO REACH THE NEAREST BOUNDARY. 


SUBROUTINE STORE 

STORES INITIAL CONDITIONS OF SECONDARY ELECTRONS FOR LATER USE BY 
SUBROUTINE HIST. IF SECONDARY STORAGE IS EXCEEDED, THE SUBROUTINE 
HALTS THE MAIN PROGRAM. 


SUBROUTINE KICK 

THE ELECTRON IS MADE TO TRAVERSE A FINAL SHORT PATH SEGMENT IN THE 
FINAL DIRECTION REACHED. THE LENGTH OF THE FINAL PATH SEGMENT IS 
TAKEN TO BE THE RESIDUAL PRACTICAL RANGE CRESIDUAL MEAN RANGE 
MULTIPLIED BY A SUITABLE DETOUR FACTOR THAT IS TAKEN FROM TRANSPORT 
CALCULATIONS OR EXPERIMENTS). 


BLOCK DATA 
DEFINES PRINTING FORMATS FOR THE OUTPUT PORTION OF ZTRAN. 


SUBROUTINE OUTPUT 
DRIVES THE OUTPUT ROUTINES OF ZTRAN. 


SUBROUTINE NORM 
NORMALIZES OUTPUT DATA. ALSO PRINTS A PORTION OF THE DATA. 


SUBROUTINE PRTOUT 
PRINTS DATA. DETAILED HEADINGS ARE PROVIDED SO THAT THE OUTPUT 
SHOULD BE SELF-EXPLANATORY. 


SUBROUTINE TAPOUT 
PROVIDES FOR WRITING OUTPUT DATA ONTO TAPE CUNIT 11). 


SUBROUTINE PUNOUT 

GENERAL NAME FOR A SUBROUTINE THAT PUNCHES OUT A DESIRED SUBSET OF 
THE DATA PRODUCED. SPECIAL VERSIONS OF THE SUBROUTINE HAVE TO BE 
WRITTEN AS THE NEED ARISES. THESE DATA CAN BE WRITTEN ONTO TAPE OR 
OTHER DATA STORAGE DEVICES BY REDEFINING THE PUNCH UNIT CUNIT 7) 
APPROPRIATELY. 
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ENERGY SPECTRA AND ANGULAR DISTRIBUTIONS OF ELECTRONS TRANSMITTED THROUGH 
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ENERGY SPECTRA AND ANGULAR DISTRIBUTIONS OF ELECTRONS REFLECTED FROM 
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